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Abstract. The problem of a particle diffusion in a fluctuating scalar field is studied. In 
contrast to most studies of advection diffusion in random fields we analyze the case where 
the particle position is also coupled to the dynamics of the field. Physical realizations of this 
problem are numerous and range from the diffusion of proteins in fluctuating membranes and the 
diffusion of localized magnetic fields in spin systems. We present exact results for the diffusion 
constant of particles diffusing in dynamical Gaussian fields in the adiabatic limit where the field 
evolution is much faster than the particle diffusion. In addition we compute the diffusion constant 
pcrturbatively, in the weak coupling limit where the interaction of the particle with the field is 
small, using a Kubo-typc relation. Finally we construct a simple toy model which can be solved 
exactly. 
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1. Introduction 

The diffusion of passive particles in complex and random velocity fields lias been extensively studied 
in statistical mechanics and fluid mechanics. In most of these studies one is interested in the 
dispersion of tracer particles which are advected by the complex or random field but which do 
not affect the field itself- it is for this reason they are called passive |l], ||. These systems can be 
studied in the case of incompressible velocity fields (relevant for problems of turbulent dispersion 
in fluid mechanics) and the case where the velocity field is derived from a gradient (relevant to 
statistical mechanics and dynamical transitions related to the glass transition). As well as diffusion 
in time dependent fields the problem of diffusion in quenched random fields has also been extensively 
studied ||. Assuming that the advecting velocity field has zero mean, the passive tracer particle 
may diffuse normally as (x 2 (£)) ~ t but in certain circumstances the particle diffuses anormally 
(x 2 (t)) ~ t 2v with v ^ 1/2. In the case where v < 1/2 the diffusion is called subdiffusive and 
when v > 1/2 the diffusion is called superdiffusive Q. Other interesting phenomena arise when one 
considers the diffusion of an ensemble of non-interacting (among themselves) tracers. Depending 
on the statistics of the advecting field clustering phenomena may arise . 

In this paper we will consider a problem where the tracer particle's position is coupled to the 
evolution of a scalar field. This means that the particle is advected by the field but also the 
dynamics of the field is affected by the particle position. We take the dynamics of the field to be 
over damped stochastic dynamics. An interesting question arises as to how the diffusion of the 
particle depends on the fields's dynamics, for instance one could have the same Hamiltonian for the 
system but in one case the field may evolve according to model A (non-conserved order parameter) 
dynamics and in the other by model B (conserved order parameter dynamics) or indeed Brownian 
hydrodynamics ||. We will first concentrate our study on the case where the field dynamics is 
much more rapid than the local diffusion of the tracer particle. In this limit we will show that the 
diffusion of the particle is always slowed down by coupling the field. This is in contrast to the case 
where the effect of the tracer on the particle is ignored and where in this limit it can be shown 
that the diffusion of the particle is speeded up. We then consider a perturbative calculation of the 
effective diffusion constant where the coupling between the field and particle position is weak. This 
computation is based on a Kubo-like relation for the effective diffusion constant, in addition this 
relation shows explicitly that the diffusion constant of the particle is reduced when its position is 
coupled to the fluctuating field. 

A concrete example of this problem is the diffusion of a protein on a membrane. The protein's 
position can be coupled to the height fluctuations of the membrane, for instance by tending to 
impose a local mean curvature Q . Alternatively the protein's position could be coupled to the local 
composition of the membrane, this could be because the protein has an affinity for a particular lipid 
type in a multicomponent lipid membrane or because it imposes a local tilt on the lipid hydrocarbon 
tails [Q . These two types of couplings lead to membrane mediated interactions between proteins 
but they will also modify the dynamics of a protein in the membrane. Indeed the question of 
what determines a proteins diffusion constant in a lipid bilayer is biologically very important. 
The problem was first addressed by Saffmann and Delbriick ^] based on a two dimensional fluid 
model. However experimental studies Jo| suggest that this simple fluid model may not explain the 
experimental data on protein and peptide diffusion constants as a function of their size. A number 
of studies have subsequently shown that protein coupling to membrane composition and height can 
substantially modify the protein's diffusion constant |l(], [yj |l2|, Q [H| [l(| . We should note that 
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height fluctuations can modify the effective observed diffusion constant of a protein in a membrane 
even when there is no coupling between the fluctuations and protein position, this is because the 
protein diffuses in the plane of the membrane but the diffusion is observed on the projected area, 
this leads to an apparent reduction of the diffusion constant, this is a geometrically induced effect 



2. The model 



Consider the dynamic of a Langevin particle whose position is denoted by x(£) diffusing with a 
linear coupling to a fluctuating free field. The overall Hamiltonian for the system is 

H =\ J 0(x)A0(x)dx - hK<f>(x(t)) (1) 

where we will take A and K to be self adjoint operators. The above Hamiltonian applies to many 
systems. For example if we take A = — V 2 +to 2 and K = 1, this is a model for a point magnetic field 
of magnitude h diffusing in a Gaussian ferromagnetic model. If A = kV 4 — erV 2 and K = —V 2 , the 
Hamiltonian is the Helfrich one for the height fluctuations of a lipid membrane where 4> represents 
the height and the choice of K is due to the fact the particle is coupled to the local membrane 
curvature . Here we are interested in the diffusion of the particle in the field. However to study 
the dynamics of diffusion we must also define the dynamics of the field. Here we will take for the 
field dynamics the general dissipative dynamics form 

and where R is a self adjoint dynamical operator and £ is a Gaussian noise of zero mean which is 
uncorrelated in time. For instance, R = 5(x — x') corresponds to model A conserved dynamics and 
R = — V 2 5(x — x') corresponds to model B dynamics [|]. If one is considering the case where <j> 
represents the height fluctuations of a membrane then using Brownian hydrodynamics the Fourier 
transform of R, R is given by -R(k) = l/47?|k|, where rj is the viscosity of the solvent surrounding 
the membrane |l7| . The field dynamics is taken so as to respect detailed balance so that the Gibbs- 
Boltzmann distribution is found for the equilibrium measure of the field and the particle position. 
This means that 

(£(x, f)£(x', t')) = 2Ti?(x - x')5(t - t'), (3) 
where T is the temperature of the system. The dynamics of the particle is given by 

-Ai = -K__ + V5*(t), (4) 

where the noise terms is Gaussian noise with zero mean zero and correlation function 

(Vi(tH(t'))=2T5 ij 5(t-t'). (5) 

The coefficients k and can be used to set the relative time scale between the dynamics of the field 
fluctuations and that of the tracer movement. In the absence of a coupling between the field and 
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the particle, the particle diffuses normally and within the notation set up here the mean squared 
displacement at large times behaves as 

(x 2 (f)) ~2dTnt = 2dDt, (6) 

where d is the spatial dimension and D = Tk is the bare diffusion constant. For the specific choice 
of Hamiltonians considered here we thus have the equations of motion 

<90(x) _ jja^) + hK(j>RK ( x _ x (i)) + y7^(x, t) (7) 



at 

and 

= kVW)) + v^fc(t), (8) 

at 

where in this paper, for two operators A and B will will denote by AB there composition as 
operators. In the limit where it is defined we will be interested in the effective diffusion constant 
for the tracer defined via 

(x 2 (i)) - 2dT Ke t = 2dD e t, (9) 

where D e is the effective late time diffusion constant. 

We note that, as mentioned above, Eq. (|J) has been extensively studied in the case where the field tj> 
evolves independently of the particle position. This problem is referred to as the advection diffusion 
of a passive scalar (the concentration of the particle) in a fluctuating field <f). This is obtained in 
the limit where one sets h — in Eq. (Q) but keeps ft, 7^ in Eq. (@). It was suggested that this 
limit can be used to approximate the diffusion of the tracer particle in |fl0| , 12 . In this case it is 



found that the effect of the field fluctuations can be to increase the diffusivity of the tracer particle 
with respect to that obtained when it is not coupled to the fluctuating field (h — 0). However 
the numerical simulations of [[l3] where the effect of the particle position on the field is taken into 
account suggests that the diffusion is reduced with respect to the case h = and in the authors 
of |l2| revisited the problem and numerically and analytically confirmed the findings of jl3| . 
The adiabatic results obtained here show that in this limit the diffusivity is always diminished with 
respect to the case h = and we argue that this is the limit where the diffusivity should be the 
most rapid. We also will show via a Kubo-type formula that active coupling to the fluctuating field 
should always reduce the value of the diffusion constant. 



3. The adiabatic limit 



We will now analyze the dynamics of this system in the limit 3> k, i.e. where the field dynamics 
is much quicker than that of the particle. The basic idea is that one can eliminate the field variable 
in a mathematically controlled manner to yield an effective diffusion equation for the particle where 
the field no longer appear implicitly. This sort of procedure can be carried out at the level of the 
Fokker Planck equation using projection operator methods |^(|. However in the case here as we 
have a dynamical variable <f> with an infinite number of degrees of freedom we will use an alternative 
method based on direct analysis of the Langevin equations |2^] . 
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We start by integrating the equation of motion of the field cj> to obtain: 

<t>(x) = ( ds exp{~ K4> (t - s)RA) [hK<f,RK(x - x(s)) + y^£(x, s)] . (10) 
Jo 

Now if «0 is large and the operator RA is positive, the above integral is dominated by the region 
where s is close to t. We make the simple change of variables u = t — s in the above to find 

</>(x) = / du exp(— n^uRA) ^hn^RKix — x(t — uj) + y^K^^x, t — u)] . (11) 
Jo 

A Taylor expansion about u = now yields for large t, 



(x)= / 



+ o(— ) 



cfej (i) 

1 

1 



exp(— k^uRA) 



h K(j ,RK(x - x(i)) + n ^ v J VjRK(x - x(t)) + ^//^(x, t) 



= ZiA" 1 ^ (x - x(i)) + —^4^{RA)- 2 V jRK{x - x(t)) + J — (i?A)^^(x, i) 
«0 at y 

(12) 

We must now compute VK(f>(x(t)) from the above. We may write the first term of Eq. ( |12| ) using 
its Fourier representation as 

hKA-'Kix - x(i)) = ^ / exp (ik • (x - x(t))) (13) 

and the second term is given by 

^^K(RA)- 2 \7 3 RK( X -x(t)) = dXl } t] £ w / rfk% - ^ (fc) exp(tk-(x-x(t))).(14) 
The results now give that to the order of approximation in l/n^ used above we have 

* R 2 (k) d Xj {t) h /■ t L k 2 (k) 



ViK<j>{x{t)) = — / dk ik^-f- £ ' \ dk kjkt^ 



(2n) d J A(k) dt k^ttY J J A{k) 2 R{k) 

+ J— V l K(RA)- 1 ^(x,t) (15) 
The first term is zero by isotropy and we can also write 

(dk k jh J^- = ^fdk eJ^B . (16) 
J J A(fc) 2 i?(fc) d J A{kfR{k) V ' 

We may thus write the effective Langevin equation for x(t) as 

Let us remark here that if we had not taken into account the effect of the particle position on the 
field and had simply considered the effect of the field on the particle we would have arrived at the 
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effective diffusion equation 



dxj(t) 
dt 



—ViK(RAy 



i £( X (t),t) + ^r !l (t) 



(18) 



The effective diffusion constant for the process of Eq. ( |l7| ) K e is simply related to that of Eq. 



1 



h 2 K 
K 4 ,d(2ir) d 



Jdk k< 



K 2 (k) 



A(k) 2 R(k) 



If we write Eq. 



as 



dxj (t) 
dt 



= G(x(t),t), 



we see that the correlation function for the noise is given by 

(C,:(x(;),t)G(x(0,t')> = 2TS(t-t')5 ijl 
which immediately yields 



K0d(27r) d 



dk k 2 -. 



k 2 {k) 



A(k) 2 R(k) 



K,pd{2TT) d 



dk k 2 ^ 



K 2 {k) 



A{k) 2 R(k) 



and 



h 2 K 
K^d(27r) 



■fdkk 2 



K 2 (k) 



A(k) 2 R(k) 



(19) 



(20) 



(21) 



(22) 



(23) 



We thus see that if the particle is not coupled to the field that the the diffusion constant (n*) is 
increased, however when the coupling is taken into account the diffusion constant is reduced. More 
over to first order in k/k^ we find that the change in the two different diffusion constants with 
respect to their bare values is the same in magnitude but of opposite sign. It is interesting to 
note that there is no temperature dependence on the renormalization of K e due to the interaction 
with the field. This means that D e retains a simple linear dependence on T within the adiabatic 
approximation. 

Note that the integrals occurring in the above can be ultra-violet (UV) divergent. This may the 
case for certain local operators K, for example _fT(x — x) = i5(x— x'). When numerically simulating 
the system in this case one can just introduce an ultra-violet cut-off in the simulation, i.e. a 
maximal Fourier mode. The result given here shows that a naive application of the Stokes Einstein 
relationship works for the effective diffusion constant D e . Define by D the diffusion constant without 
a coupling to the field. For a particle moving at constant velocity the frictional force, opposing the 
motion, is given by 



/o = A w 

However, the friction is related to the mobility via 

Mo = fov 



(24) 



(25) 
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Stokes Einstein tells us that, when it is valid (see later discussion), 

D = MoT, (26) 
and using the fact that D = Tk we have that Ao = 1/k 

Now in presence of the coupling one can compute the average frictional force due to the 
fluctuation field [|l4| [l5| and one finds that 

U = V ( 27 ) 

where A^ is given by 

A , = - / dk k 2 . R2{k) (28) 

K(j> d{27r) d J A(fc) 2 i?(/fc) V ' 

The total frictional force is thus given by 

/ = fa + U = (Ao + A*)u, (29) 

This gives via the Stokes-Einstein relation 

T 

D e = fi e T = 

*a + *<t> 

D 

1 I h?Df} \ frfk h 2 R2{k) 1 ' 

where (3 = 1/T, which is equivalent to the result Eq. (|23). We also note that all the terms in 
the integrand of Eq. (|2^) are positive and thus we have within the adiabatic approximation that 
D e < D. The application of the Stokes-Einstein relation that we have just made is clearly not exact. 
To compute the diffusion constant using the Stokes-Einstein relation one must compute the average 
value of the velocity v at constant applied force || . This is a much harder problem than computing 
the average force at constant velocity. In a previous paper [[l5| we argued that the Stokes Einstein 
relation as applied above should be valid when the fluctuations of the force are small and thus the 
force is near to constant in the statistical sense. The fact that the average force is large means that 
the friction is large and thus the diffusion constant is small. Here we see that the adiabatic limit 
reproduces the approximate application of the Stokes-Einstein relation given above. This result 
can be explained by examining the expression given for the autocorrelation function of the force 
fluctuations given in |l4| , here one sees in the adiabatic limit that fluctuations of the force become 
uncorrelated in time and their amplitude becomes small. 

3.1. Examples 

We will now consider a number of special cases of our principle result Eq. (^3|) which we will write 
as 

K 



i , h 2 S d KQ 
1 K ^d(2Tr)" 



(31) 
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with 



Q= dk fc a+i -= M— , (32) 

and where Sd is the area of a sphere of radius 1 in d dimensions. 

Considering now a localized magnetic field diffusing in a ferromagnet within the Gaussian 
approximation, we take 

A(k) = k 2 +m 2 (33) 

we assume a magnetic field with a localized Gaussian profile and thus 

/ h 2 n 2 \ 

K(k)=exp{-—). (34) 
With this we find 

roc d+1 

Q = m d - 2 J o d 1 {q l +l) 2 cxp(-g 2 m 2 a 2 ) (35) 



for model A dynamics and 

/>00 

/ df, 

(q 2 + l) 



Q = m d - 4 dq f I ; exp(-g 2 m 2 a 2 ) (36) 
Jo 



for model B dynamics. 

In certain cases one can take the limit a — > in the above and thus obtain results that are only 
weakly dependent on the cut-off. However this cases depend strongly on the mass of the field theory 
and the result can be divergent when the theory is critical, i.e. when m = 0. These cases are the 
following 



Model A: d < 1 



7T 

4m' 



d = 1. (37) 



Mode/ 5: d < 3 



Q 4m 3 ' 



f 



d = 2 

d = 3. (38) 



2m 2 ' 
~ 4m' 

We thus see that as the mass of the scalar field is decreased or its correlation length £ = 1/m 
increases the diffusion constant of the active tracer particle is decreased. In the limit where 
fi > 1 the magnetic tracer will be surrounded by a polarized region where the field <f> has the 
same sign as the tracer field. The size of the polarized region will be of order £ and the modification 
to the diffusion constant above presumably reflects the effective mobility of this polarization 
cloud. It is also interesting to note the different dependence on m between model A in B in 
one dimension. 
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Cases where the results have a strong dependence on the cut-off have been extensively discussed in 

4. Weak coupling limit 

In this section we use a formally Kubo-like expression for the effective diffusion constant, the 
expression can formally be computed to 0(h 2 ) in the particle-field coupling parameter, thus giving 
an expression for D e which is exact to this order. To obtain the Kubo formula we integrate Eq. (j|) 
to obtain 

x(t) - x(0) = hK [ VA>(x(s))ds + VWnB(t) (39) 
Jo 

where B t is a standard d dimensional Brownian motion with 

(B t {t)B 3 {s)) = 5ij min{t, s}. (40) 

In the above we assume that at t = the system is in equilibrium (i.e. we assume that the volume 
of the system is finite and Xo and the initial field configuration </> are chosen from the equilibrium 
distribution (see || for more details) . Now subtracting the first term of the right hand side of Eq. 
( p9| ) from both sides, squaring the resulting equation and taking the average yields 

((x(i) - x(0)) 2 ) + h 2 n 2 / ^ VK(f>(x(s))ds^ \ = 2Tdnt, (41) 

where the cross term on the right had side is zero due to the Onsager relations ||. We can thus 
define a time dependent diffusion constant via 

((x(i)-x(O)) 2 ) = 2dD e {t)t = 2dTn e (t) = 2Tdnt-h 2 / J ^ VK4>(x(s))ds^ \ (42) 

where the late time limit of these two quantities are the effective values limt-yoo D e (t), K e (t) = 
D e , K e . We may therefore write 

t^-^fi'i(Ki' VA * w)< ' s ) 3 )^ (43) 

^From this exact formula we see that the value of the diffusion constant is reduced by the interaction 
with the field. This expression may now be evaluated to 0(h 2 ) by replacing x(t) by the pure 
Brownian motion \/2-DB t of the particle without interaction of the field and using the correlation 
function for the free field without interaction with the particle which can be written as 

(0o(x, t)to(y, s)) = T J J^-A-^k) exp (-«^|t - s\R(k)A(k)) exp (ik • (x - y)) .(44) 

After a straight forward computation using the fact that B(i) and 4>q are independent we 
obtain 

D '-l-!*>°[ A,-- , ekHk) r + 0(h% (45) 



D 



( 2n ) d A(fc) (^A(fc)A(fc) + Dk 2 
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Figure 1. Effective diffusion coefficient versus field evolution velocity for the model of Eqs. ( 
and (ttfl) in one dimension for V(x) = cos(a;) and temperatures T = 0.5 (+), 1 (x) and 2 
The continuous lines correspond to the analytical result. 



Note that in the adiabatic limit this result is clearly equivalent to Eq. ( J30| ) to 0(h 2 ). In terms of 
the variables re Eq. (|4^) reads 

h 2 n f dk k 2 K 2 (k) , , 

a .1 [Zir) a a fu\ I »■ . vth\ a tv\ _i_ Tnk 2 



( 2n ) d A(jfe) ( K4> R{k)A(k) 



note therefore in contrast to the purely adiabatic result Eq. ( |23| ) that there is a temperature 
dependent renormalization of re e . Another interesting thing to note is that in this perturbative 
result we can recover the case where the field is frozen, i.e. the other extreme to the adiabatic limit 
where re^ = 0. Here the field is quenched and has correlation function TA _1 . This quenched result 
agrees with the first order perturbation result for quenched random fields §|. 



5. A Toy Model 

In the general class of models studied above we have been able to obtain partial results on the 
effective diffusion constant on an active tracer in two distinct limits, the adiabatic limit and the weak 
coupling limit. Here we present a simple toy model whose behavior can be thoroughly analyzed. 
We consider a tracer particle x coupled to another diffusion process y via the Hamiltonian defined 

by 



ff(x,y) = V(x-y), 



(47) 
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where V is a function such that there exists a vector (or set of vectors) a such that, V(x) = V(x + a) 
in the algebraic or statistical sense. For example we could take V(x) to be a periodic function or 
one that is statistically invariant by translation. The coupled diffusion equations for x and y are 
given by 

dxi dV(x — y) , . . 

- ' +V^Vx t (48) 



dt dxi 
dyi dV(x — y) 

~ul or, ; ^~' / "- (49) 

where the noise variables above are white noise at temperature T as defined earlier. In the case where 
the variable y is frozen (or equivalently k v = 0), as the function V is translationally invariant and 
if it is bounded, we expect the process x will have an effective quenched diffusion constant defined 
by 

(at?) ~ 2dD{ q) t = 2dTn ( e q) t, (50) 

and which is independent of the choice of y by the translational invariance of V. We note that it 
is possible to compute D e exactly in a number of special cases || . In order to see the effect of the 
dynamics of y on the process x we define the new variables 

u = x-y (51) 

v = K y x + n x y, (52) 
and it is then easy to see that these new dynamical variable obey 
du t _ ~.9V(u) — 

~~q£ — — \ Kx ' K v ) — Q~a K xflxi — v K yVyi (5o) 



~~q£ = K y\/ K xVxi + K x\/ K yVyi ■ (54) 
Furthermore one can easily see that the noise in these two equations are independent and thus the 
processes u and v are independent. The process u is simply a time rescaled form of the quenched 
problem and v is a free Brownian motion. The mean squared displacement of the two processes 
can thus be computed easily and are given by 

(u?)~2dT/4«>^±^i (55) 

(vj) ~ 2dT k x k v (k x + K y )t. (56) 
Finally using the independence of u and v we find that 

(x?) ~ 2dTn x f'^J^y) t , (57) 

\ K x + K y I 

which gives the effective diffusion constant of x as 

D e ( + Ky\ , . 

-f = ^ = "A K +K • (58) 

J V l^x T riy I 

We see that in the quenched limit n y — we obtain that n e — as we should. In addition 
it is clear that n e is an increasing function of K y , the quenched case being a lower bound for the 
effective diffusion constant. Another interesting fact about the expression Eq. (p8[ ) is that there are 
a number of models where the quenched diffusion constant vanishes signaling the transition from a 
regime of normal diffusion to one of subdiffusion [|, ||^, Q , however we see from Eq. (|58|) that 
when K y 7^ then the vanishing of does not cause the diffusion constant to vanish as the result 
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has an additive property. Indeed if «4 = then we find 



U i I I V I J \ V V ■ I 1 L I ' I 

(59) 

Ky X ~\- Ky 

throughout the parameter region where the quenched problem shows subdiffusion. 

The quenched diffusion constant can be computed exactly in a limited number of cases [0, 
notably in one dimension where it is given by 



& 1 



k x (exp(-0V))(exp(pV)) 

where 



(exp(±pV)) = Km ~ [ exp(±(3V{x))dx 



(60) 



(61) 



which exists for translationally invariant potentials. In addition when the field V is equivalent to 
— V (either functionally or statistically) a exact result in two dimensions gives 



1 



k x (exp(l3V)) ' 



(62) 



As a test of this result in one dimension we took the potential V(x) — cos(x) and Eqs. (58) and 
( |60[ ) then give 

( ifm + Ky \ 



- , (63) 

\ K x "T Ky I 

where Iq(z) denotes the modified Bessel function. Numerical simulations were carried out for k x = 1 
and K y varying between 0.1 and 10 and at temperatures T = 0.5, 1 and 2. The results given Figs. 
(Q) show an excellent agreement with the analytical result Eq. (|3|). 



6. Discussion 



We have examined the diffusion of an active tracer particle coupled to a fluctuating field. Most 
previous studies have been carried out on passive tracers in time dependent or quenched fields. The 
action of the tracer on the field means that diffusion is always slowed down with respect to the non- 
interacting case, this fact is explicit in Eq. (|43|). This result may seem odd from a physical point 
of view as one would naively expect that the fluctuating field would help the particle to diffuse. 
Indeed we have seen that in the adiabatic approximation a passive diffuser diffuses more quickly 
when driven by the field. However when the effect of the tracer on the field is taken into account the 
effect of the extra noise from the fluctuating field is eliminated from a additional drag due to the 
action of the tracer on the field. Although it is in a limiting case the adiabatic calculation carried 
out here shows how these two effects compete and lead to slowing down of the diffusion with respect 
to the free case. As an example we analysed the diffusion constant of magnetic fields diffusing in 
Gaussian ferromagnets, here we found that the diffusion rate decreases as the correlation length of 
the field increases. This result is presumably linked to the fact that the field is trapped in domains 
where the field has the same sign and also its presence leads to the formation of these domains 
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about it. As the correlation length increases the size of the domains containing the held increases 
and the diffusion constant of the held effectively becomes that of it surrounding domain. It seems 
physically reasonable that the diffusion constant of the surrounding domain becomes smaller as its 
size increases. 

We have also analysed the active diffusion problem in the weak coupling limit where the existence of 
an underlying Gibbs measure allows us to write down a Kubo-type formula for the diffusion constant. 
As well as rigorously establishing that diffusion is slowed down with respect to the free case, this 
formula can be used to give a hrst order expression for the modified diffusion constant. 

Finally we have analysed a toy model for a particle interacting with a scalar held in the simple 
case where this scalar held is another diffusing particle. In this case the effective diffusion constant 
can be formally computed in terms of the effective diffusion constant for a particle diffusing in a 
quenched random potential. The effect of the dynamics of the second diffuser can be thoroughly 
understood and as its bare diffusion constant is increased so is that of the tracer. 

Clearly there are still a large number of questions about diffusion of active scalars in fluctuating 
holds, for the precise problem examined here the whole regime beyond the perturbative and 
adiabatic regimes studied here remain open. It would also be interesting to understand in more 
detail the cross over from the active to passive cases and understand under what circumstances 
held fluctuations increase/decrease the tracer diffusion constant. 
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